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The communication capacity of Gaussian bosonic channels with memory has recently attracted 
much interest. Here, we investigate a method to prepare the multimode entangled input symbol 
states for encoding classical information into these channels. In particular, we study the usefulness 
of a Gaussian matrix-product state (GMPS) as an input symbol state, which can be sequentially 
generated although it remains heavily entangled for an arbitrary number of modes. We show that the 
GMPS can achieve more than 99.9% of the Gaussian capacity for Gaussian bosonic memory channels 
with a Markovian or non-Markovian correlated noise model in a large range of noise correlation 
strengths. Furthermore, we present a noise class for which the GMPS is the exact optimal input 
symbol state of the corresponding channel. Since GMPS are ground states of particular quadratic 
Hamiltonians, our results suggest a possible link between the theory of quantum communication 
channels and quantum many-body physics. 
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I. INTRODUCTION 

Quantum communication channels are at the heart of 
quantum information theory. Among quantum channels, 
the bosonic Gaussian channels describe very common 
physical links, such as the transmission via free space or 
optical fibers pQ . The fundamental feature of a quantum 
channel is its capacity, which is the maximal information 
transmission rate for a given available energy. The capac- 
ity can be classical or quantum, depending on whether 
one sends classical or quantum information (here, we fo- 
cus on the former). In previous works, it was shown 
that for certain quantum memory channels, in particular 
channels with correlated noise, the optimal input symbol 
state is entangled across successive uses of the channel; 
see Refs. [2TfTT] and references therein. In general, such 
multimode entangled states may be quite hard to pre- 
pare, which motivates the present work. 

In this paper, we address the problem of implementing 
the (optimal) input symbol state for Gaussian bosonic 
channels with particular memory models. For this pur- 
pose, we study the usefulness of the so-called Gaussian 
matrix-product state (GMPS) [TH H3] as an input sym- 
bol state for the Gaussian bosonic channel with addi- 
tive noise [HI HJJ and the lossy Gaussian bosonic channel 
[3 HI [10]. This translationary-invariant state is heav- 
ily entangled and can be generated sequentially, which 
happens to be crucial for its use as a multimode input 
symbol state in the transmission via a Gaussian bosonic 
channel. The GMPS are known to be a useful resource 
for quantum teleportation protocols [THUS], but, to our 
knowledge, they have never been considered in the con- 
text of quantum channels. 

In Sec. II, we give an overview of the method used to 
derive the Gaussian capacity of Gaussian bosonic mem- 
ory channels, following our previous work [3][S][TT]. Our 
original results are presented in Sec. Ill, where we ad- 
dress the use of GMPS in this context. In Sec. Ill A, we 
show that the GMPS, though not being the optimal input 
state, is close-to-capacity achieving for Gaussian bosonic 



channels with a Markovian and non-Markovian noise in a 
large region of noise correlation strengths. In Sec. Ill B, 
we provide a class of noisy channels for which the GMPS 
is the exact optimal input state. Since the GMPS is as 
well the ground state of particular quadratic Hamiltoni- 
ans, this suggests a direct link between the maximization 
of information transmission in quantum channels and the 
energy minimization of quantum many-body systems. In 
Sec. Ill C, we also observe that the squeezing strengths 
that are needed to realize the GMPS in an optical setup 
are experimentally feasible. Finally, our conclusions are 
provided in Sec. IV. 



II. GAUSSIAN CAPACITY OF MEMORY 
CHANNELS WITH CORRELATED NOISE 

A. Gaussian bosonic channels 

Let us now consider an n-mode optical channel fW, 
which can either be a bosonic additive noise channel or 
a lossy bosonic channel. In the following, n single-mode 
channel uses will be equivalent to one use of an n-mode 
parallel channel PJj] . Each mode j is associated with the 
annihilation hj and creation operators, or equivalently 
with the pair of quadrature operators qj = (&j + <zt)/-y/2 
and Pj = — a,j)/y/2, which obey the canonical com- 
mutation relation [qt,Pj] — i5ij. By defining the vector 
of quadratures R = (q%, q n ;pi, ■ ■■,Pn) T , we can ex- 
press the displacement vector m = Tr[pH] of any state 
p, along with its covariance matrix (CM) 



with 



7 = Tr[(R -m)p{R- 



m) T ] - J/2, 



where I is the n x n identity matrix. In phase space, a 
Gaussian state is defined as a state p having a Wigner 
distribution that is Gaussian; hence, it is fully character- 
ized by its mean m and CM 7. 
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For the channel encoding, we consider a continuous 
alphabet, that is, we encode a complex number q + ip 
instead of a discrete index into each symbol state. We 
encode a message of length n into a 2n-dimensional real 
vector r = q 2 , ... , q n ;Pi,P2, — ,Pn) 7 - Physically, this 
encoding corresponds in phase space to a displacement 
by r of the n-partite Gaussian input state defined by its 
mean mj n and CM 7; n . The modulation of the multipar- 
tite input state is taken as a (classical) Gaussian multi- 
partite probability density p mo d( r ) with mean m mo( j and 
CM 7 mo d- The means of the input state wij n and clas- 
sical modulation m mo( j can be set to zero without loss 
of generality because displacements leave the entropy in- 
variant; hence, they do not play any role in the capacity 
formulas defined in Sec. II B. The action of the chan- 
nel T*-™' is thus fully characterized in terms of covariance 
matrices, that is, 

Tout — ^ Tin ^ 7env 5 
7 Tout ~^ ^ 7mod ; 

where 7 ou t and 7 are the CM of the individual out- 
put and modulated output states, respectively. For 
K = k' = 1, Eq. ([l]) defines the bosonic Gaussian channel 
with additive noise, where 7 cnv is the CM of a (clas- 
sical) Gaussian multipartite probability density p ell v( r ) 
describing noise-induced displacements in phase space 
(see Ref. [TT] for details). For k = r\ and k? = 1 — n, 
with a beamsplitter transmittance r\ € [0, 1], Eq. ([!]) de- 
fines the lossy channel where 7 env stands for the CM of 
the environment state (see Ref. [TO] for details). Both 
channels obey the physical energy constraint that reads 
Tr(7i n + 7mod)/(2n) — 1/2 = n, where n is the mean 
photon number at the input. 

B. Gaussian capacity 

In recent works, we found the Gaussian capacity (i.e., 
the capacity when restricted to Gaussian input states ac- 
cording to the usual Gaussian channel minimum entropy 
conjecture) and optimal input encoding for the additive 
Gaussian channel with noise correlations between subse- 
quent uses of the channel modeled by the CM [5J [TT] 

7cnv = ( 7c o llv iL) (2) 

where 7| nv and 7 p nv are commuting matrices of dimen- 
sion n x n. The absence of correlations between q and 
p in Eq. ^ is generally considered to describe a natural 
noise. We found that the optimal input and modulation 
CM 7 ; * n and 7mod are diagonal in the same basis as the 
noise CM 7 onv , and have the same block structure. Thus, 

7f„ = 7ta © ^in and 7mod = 7mod © T&od- In addition, 
the optimal input state is pure, i.e., det (27;* n ) = 1, which 
implies 



From now on, we consider the optimal input and modula- 
tion eigenvalue spectra in the limit of an infinite number 
of channel uses n — > 00, so all matrices must be expanded 
to infinite dimensions, see Ref. [5]. 

For an input energy n above a certain threshold 
n t hr, the optimal eigenvalue spectra are linked via 
a global quantum water filling solution [9], that is, 
7 <3 *(a;) = 7 p *(x) = const., Va; € A where x is a contin- 
uous spectral parameter within a spectral domain A and 
7 ?,p * (x) is the spectrum of the q and p blocks of the opti- 
mal modulated output CM 7* = 7 9 *©7 P *. Furthermore, 
the optimal input state was determined as [HI fTU] 

7m { } 2y 7 iw [ ) 

which corresponds more precisely to the spectrum of the 
q and p blocks of the optimal input CM, 7* n . We remark 
that this holds for both the additive noise [H] and lossy 
channel |10) . 

In the following, we will consider noise models (see Sec. 
II C) characterized by a CM with symmetric spectrum, 
i.e., 7| nv (a;) = 7£ nv (|-4| — x), where |^4| is the size of 
the spectral domain A. Furthermore, the noise models 
fulfill maxa;{7| nv (a;)} = 7f nv (0). For this case, the input 
energy n that is required to fulfill the global quantum 
water filling solution and Eq. Q, for all x, is given by 

n > nt„ r = it (0) + 7lnv (0) - 5 - (5) 

where N — A J xeA dx 7onv( a; ) stands for the added noise 
energy. Throughout this paper, we only consider the case 
above threshold, when n > n tn r- Then, the Gaussian 
capacity of the channel with additive noise is given by 

m 

C = g(n + N) 

~R| / dxg (v 7 ™ 4 ^ 7 ™^ ~ 5)' ( 6 ) 

xeA 

where 7out*( a; ) = 7i„ P *( a; ) +7onv( ;E ) according to Eq. ([j}. 
The function g(x) stands for the entropy of a ther- 
mal state with x photons. It is defined as g(x) = 
(x + 1) log (x + 1) — x log x if x > 0, and g(x) = if 
x < 0, where log(a;) denotes the logarithm to base 2. 

Now, if one restricts the input states to independent 
coherent states in the case of global water filling, one 
may also define the coherent- state rate which is 

given by Eq. ^ replacing 7^* (x) by 1/2 + 7^(2;). 
For the lossy channel the quantities given by Eqs. ^ 
and ^ as well as the coherent-state rate are obtained 
by replacing n rjn,^ v {x) -4- (1 - vhcn V ( x ) and 
7in P *( a ~* ? 77in P *( x )- Note that all these expressions 
also rely on the assumption that the Gaussian capacity 
of independent Gaussian channels is additive, see Ref. 
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Noise models 



(a) 



Let us introduce two different noise models which will 
be used to model the Gaussian memory channels, namely 
a Markovian and non-Markovian model. 



1. Markov additive noise 

In Refs. [HUH], we considered a classical Markov noise 
with variance Nm > 0, given by 



7envM = Nm 



M{4>) 








(7) 



where M(<j>) is an 



n x n 



matrix defined as Mij((f>) 



^h-il^ w ith the correlation parameter Q < <f> < 1. Note 
that M((j>) and M(— <fi) commute in the limit of an infi- 
nite number of channel uses. In this limit, the spectra of 
the quadrature blocks 7eAvM = NmM(±4>) are given by 



i> 2 =F 2ipcos(x) 



arG[0,27r], (8) 



with the upper (lower) sign standing for the q (p) quadra- 
ture. By using Eq. Q, we find that the optimal input 
state is an infinite product of squeezed states. Then, 
when rotated back to its original basis, the optimal in- 
put state becomes a multimode entangled state [9l E] . 



2. Non-Markovian noise 

A non-Markovian channel noise model was considered 
in Refs. [1(10], given by 



7onvAT - N N I q 



,sf2 





-sfl I J 



(9) 



with iVjy > 1/2 for the considered lossy channel (N^ > 
for a non-Markovian additive noise channel), s € K, and 



where 1~2 is a n x n matrix defined as fL 



V~i.j + l+0~i+l,j ■ 

N N e ±sn 



The spectra of the quadrature blocks "YciivN 
read 



^ vN (x)=N N e ±2sc °<*\ ie [0,H (10) 

with the upper (lower) sign standing for the q (p) quadra- 
ture. In the case of a global water filling, it was shown 
that the optimal input state [Eq. Q] is also entangled 
in the original basis [51 [TU] , as for the Markov additive 
noise. 

Since the optimal input state for both noise models 
exhibits multimode entanglement across the subsequent 
uses of the channel, with n — > oo, its preparation may be 
a very challenging task. This is what we investigate in 
the next section. 
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FIG. 1. (a) Optical scheme for the preparation of the Gaus- 
sian matrix-product state (GMPS), slightly modified with re- 
spect to Ref. [13] . Here, TMS stands for a two-mode squeezed 
vacuum state with squeezing r-r, while 7s represents the 
three-mode building block. Note that all TMS and three- 
mode building blocks could each be generated by a single de- 
vice that is used repeatedly. One half of the TMS generated 
at time i is used immediately to generate the GMPS mode 
i, while the other half is sent to a delay line (to be used at 
time i + 1). After two Bell measurements (represented by 
curly brackets) involving the two TMS halves (noted I and r) 
and the two upper modes of fB (noted 1 and 1') followed by 
appropriate conditional displacements, the third mode (noted 
2) of ~fs collapses into the GMPS mode i. (b) Optical setup 
of the three-mode building block -ys that is used to generate 
a nearest-neighbor correlated GMPS. Here |0) denote vacuum 
modes, S(tb) is a one-mode squeezer with parameter r_g, and 
the bold horizontal bars represent 50:50 beamsplitters. 



III. GAUSSIAN MATRIX PRODUCT STATE 

We now address the question of how to optically imple- 
ment the optimal input states. In this context, we exam- 
ine the so-called Gaussian matrix-product state (GMPS), 
which is heavily entangled just as the optimal input state, 
has a known optical implementation, and can be gener- 
ated sequentially. This state was first discussed in Ref. 
[12"] as the ground state of particular Hamiltonians of 
harmonic lattices. In general, GMPS are constructed 
by taking a fixed number M. of finitely or infinitely en- 
tangled two-mode squeezed vacuum states shared by ad- 
jacent sites, and applying an arbitrary 2M. to f mode 
Gaussian operation on each site i. 

In what follows, we restrict our discussion to a pure, 
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translationally-invariant, one-dimensional GMPS, and, 
furthermore, to a single finitely entangled two-mode 
squeezed (TMS) vacuum state per bond between adja- 
cent sites (A4 = 1). We use the protocol introduced in 
Ref . [13] , depicted in a slightly modified form in Fig. [lja) . 
Each GMPS mode i is obtained by operating on a three- 
mode entangled state (called "building block" ; see Refs. 
[TBI ITS] for details) together with the shares (/ and r) 
of the two TMS vacuum states connecting site i to the 
left and right sites, respectively. As shown in Fig. [lja), 
a first tcleportation is performed by making a Bell mea- 
surement on modes I and 1, followed by a conditional 
displacement on mode 1'. A second teleportation then 
is made with a Bell measurement on modes r and 1', 
followed by a conditional displacement on mode 2. The 
final state of mode 2 then reduces precisely to that of 
the ith mode of the desired GMPS. We focus now on the 
mathematical description of the GMPS and its use as an 
input state, while we discuss its experimental realization 




with single- mode squeezers in Sec. Ill C 



FIG. 2. Rates of a channel with additive Markov noise: Gaus- 
sian capacity C (solid line), GMPS-rate Rgmps (crosses) and 
coherent-state rate R co h (dashed line) vs. correlation <j>, where 
from top to bottom Nm = {0.5, 0.7, 1}. We took n — 5. 



A. GMPS as approximating input state 



The CM of the GMPS can be written as 



f GMPS — 



i (c 







(11) 



where C is a n x n circulant symmetric matrix. In 
Rcf. [12 , it was proven that the correlations of a one- 
dimensional GMPS decay exponentially. Therefore, in 
the limit n — > oo, the spectrum of C _1 reduces (up to a 
change of variance) to the spectrum of M(<ft) [5T], that 



1 



A (c 1] (x) = 7gmps( 2: ) 



N 



1 + 0in - 2 4>in COs(x) 



+ A ,(12) 



with x G [0,2tt], N > 0, < <p in < 1, A G M, and 
the additional condition AN > —1/2 ensuring that the 
spectrum corresponds to a quantum state. By com- 
paring the spectrum of Eq. (12) with the optimal in- 
put spectra [Eq. (4])] for the noise models of Eqs. Q 
and one can directly verify that the optimal input 
state is not a GMPS. However, one may use the GMPS 
as an approximation of the optimal input state for both 
these noise models. By calculating the transmission rates 
for noise models [Eqs. and ^} with the GMPS 



as input state [using Eq. (p 
7 G mps( x ) + 7&v(^)]) we nnc 



and replacing ^*{x) by 
numerically that the high- 



est transmission rate is achieved for a GMPS with near- 
est neighbor correlations 7gmps n.n. [H] • We find that 
among all GMPS given by Eq. which can be gen- 

erated with the setup defined in Fig. [T] only the GMPS 
with nearest neighbor correlations has a symmetric spec- 
trum, that is 



7GMPS,n.n.( a; ) — ^GMPS.n.n.l 71 " x )- 



(13) 




FIG. 3. Rates of a channel with additive non-Markovian 
noise: Gaussian capacity C (solid line), GMPS-rate -Rgmps 
(crosses) and coherent-state rate R co h (dashed line) vs. cor- 
relation s, where from top to bottom iVjv = {0.5, 0.7, 1}. We 
took n = 5. 



Since the noise spectra defined in Sec. II C satisfy the 



same symmetry, it is intuitively clear that this type of 
GMPS is the most suitable state for these noise mod- 
els. The optical setup for the three-mode building block 
that generates this nearest-neighbor GMPS is depicted 
in Fig. [IJb). More details on it are provided in Sec. Ill 
C. 

From Eq. |l3| and the fact that the GMPS used as an 
input is a pure state, i.e., 7gmps(- c ) 7gmps( 2 ') = 1/4, Vx, 
we find that N = (1 + </f n )/(l - 0f n ) and TV A = -1/2. 
Thus, the nearest neighbor correlated GMPS has quadra- 
ture spectra 



Tgmps ,n.n. 

(x) 



1 



1 



(14) 



1 + tfn T ^^in COs(x) 2' 

with the upper (lower) sign for the q (p) quadrature. 
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FIG. 4. Rates of a lossy channel with non-Markovian noise: 
Gaussian capacity C (solid line), GMPS-rate -Rgmps (crosses) 
and coherent-state rate R co h (dashed line) vs. correlation 
s, where from top to bottom r\ = {0.5,0.7,0.9}. We took 
Nm = 1 and n = 5. 
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FIG. 5. Optimal input correlation <j}* n (solid line, left axis) 
and corresponding squeezing n n (dashed line, right axis) vs. 
correlation {<f> or s) for (a) the channel with additive Markov 
noise, where the crosses depict 0/2; (b) the channel with non- 
Markovian noise (lossy and additive), where the crosses depict 
s/2. We took for both plots Nm = Nn = 1 and n — 5. 



Therefore, when looking for the optimal transmission 
rate, one has to optimize only over the parameter (f>- m . 
In order to satisfy the global water filling solution for 
the GMPS, we replace 7i 9 n *(0) by 7 G MPS,n.n.(0) m Eq. (§, 
which leads to a modified input energy threshold depend- 
ing on </>;„, that is, 



-GMPS 
'thr 



nt hr -[7? n *(0)-7 G MPS,„.n.(0)]- (15) 



As we require that the input energy n > n^ PS , Eq. ( 15 ) 



imposes an upper bound on (f>- ln . 

In Figs. [2]|4j we plot the rates obtained for the GMPS 
with the spectrum given by Eq. ( 14 1 calculated via a 
maximization over 4>- ln , which we denote as -Rgmps- I n 
Fig. [2] we observe that for the channel with additive 
Markov noise -Rgmps is close-to-capacity achieving; 
in the plotted region, -Rgmps/C > 0.999. For the addi- 
tive channel with non-Markovian noise ([9]), we conclude 
from Fig. [3] that the GMPS serves as a very good re- 
source as well; in the plotted region, -Rgmps/C > 0.999. 



We confirm the same behavior for the lossy channel with 
non-Markovian noise, as shown in Fig. [4] for different 
beamsplitter transmittances 77. 

The optimal input correlations </>* n for both noise mod- 
els are approximately given by <p/2 and s/2, respectively, 
as can be seen in Fig. ga) and Fig. [gb). This can be 
verified as follows. Since the quantum water filling solu- 
tion holds for the GMPS with nearest neighbor correla- 
tions, its rate is given by Eq. ^ replacing Jou t *( x ) by 
7gmps n.n. ( x ) + linvi x )- I n order to find the optimal <rj in 
it is sufficient to minimize only the second term in Eq. ^ 
as only this term depends on <j)- m . This term is a definite 
integral of a function whose primitive is not expressed in 
terms of elementary functions and <pi n . However, if the 
integrand as a function of parameter </j; n can be properly 
minimized for all values of the variable of integration x 
the integral will also be minimized. In order to verify this 
possibility we take the first derivative of the integrand 
and set it to zero. This leads to the following relation: 



7e V(aO (l + 0*„ 2 + 2Ccosa;) 2 



7e P nv(z) (l + 0* n 2 ~20* n COSx) 2 



(16) 



As it happens in the general case, there is no unique pa- 
rameter c/>* n which satisfies Eq. (16) for all x. Neverthe- 
less, it is possible to obtain an approximating equality 
by neglecting the quadratic and higher order terms in 
the noise spectra given by Eqs. §7§ and ^ and in the 
right-hand side of Eq. (fl6|), i.e. 



1 + 2acos(a;) 1 
1 — 2a cos(a;) 1 



m n cosjx) 

40* n cos(a;) 



(17) 



where a = 4> for the Markovian noise and a — s for 
the non-Markovian noise, respectively. This is a valid 
approximation taking into account that cf>* n < 1 and can 
be satisfied by a unique parameter (j)* ln for all x. Namely, 
we find the simple relations cj>? ~ <fi/2 and cf>* n s» s/2, as 
verified in Fig.[5]ja) and Fig.JHJb), respectively. 



B. GMPS as exact optimal input state 

Although we have seen that the GMPS is not the op- 
timal input state for the noise models introduced in Sec. 
II C[ it is possible to do better. Indeed, for all noises 
given by 



7cnv = (A/"env © A/" nv) X (C 1 © C), 



(18) 



where Af onv is an n x n matrix that commutes with C 



given in Eq. (Ill, the GMPS is the exact optimal input 
state, that is 



7*n = 7GMPS, n>n 



GMPS 
thr ' 



(19) 



where now trivially n£j^ PS = n t hr- This is a direct result 
that can be deduced from the shape of the CM 7 gmps 
and the fact that the CM of the optimal input state (given 



G 



by Eqs. ^ and Q) is diagonalized in the same basis as 
the CM of the noise. 

Furthermore, as already mentioned, GMPS are known 
to be ground states of particular quadratic Hamiltonians 
[12] . More precisely, 7gmps is the CM of the ground 
state of the translationary invariant Hamiltonian, given 
in natural units by 



H 



(20) 



where qi and fn are the position and momentum opera- 
tors of an harmonic oscillator at site i and the potential 
matrix is simply given by V — C 2 , where C is defined in 
Eq. pi. 



A realistic example for a noise of the shape of Eq. ( 18 ) 
is given by the CM of the (Gaussian) state of the system 
defined in Eq. (201, i.e., a chain of coupled harmonic os- 



cillators at finite temperature T. We assume the system 
to be described by a canonical ensemble, thus the density 
matrix of the oscillators is given by the Gibbs-state 



PG = 



exp {-PH) 
Tr[cxp {-PH)] ' 



(21) 



where ft = 1/T. The CM of the Gaussian state pa is 
given by Eq. ([l8| withA/" cn v = I+[2 exp (/3C) — I]" 1 (see 
Ref. [19] for details), where indeed [Af cnv ,C] — 0. There- 
fore, if we assume the noise of the channel to result from 
a chain of coupled harmonic oscillators at finite temper- 
ature T, that is, 7env = 7g, then the GMPS with CM 
7gmps is both the ground state of the system given by 
Eq. (20) and the exact optimal input state for n > n t h r - 



C. Experimental realization 

Let us finally discuss the required optical squeezing 
strength to realize the optimal input correlation <p* n for 
both noise models. We first present the mathematical 
description of the three-mode building block that gener- 
ates the GMPS with nearest neighbor correlations. The 
CM of this building block is given by [TB] 
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(22) 



with w = (i + l)/2, v = (t - l)/2 and u = ^(t 2 - l)/2, 
where t > 1. The optical scheme for the three mode 
building block is depicted in Fig. [TJb) , where S(rs) is 
a one-mode squeezer with parameter tb such that t = 
cosh(2rs) [T5]. The resulting CM of the n-mode pure 
GMPS is given by [T3] 



7GMPS 



6T 



TMS C 



lp 



(23) 



with 



-I, where / is the n x n identity matrix 



1 n 

r * = 2® dia s{*>*}' 

r T _ 1 /T\ fu u 
1 wt - 9 tf 



2 vl/ v -u -u 

i=l 



(24) 



1 2n / 
(_T_i ( w 



w v 
w 



TtMS = ^TMsifT) © 7TMs( — Tt), 



where 7tms(»*) 

/ch(r) 

ch(r) sh(r) ••• 

sh(r) ch(r) ••• 

: ch(r) sh(r) 

: sh(r) ch(r) 



sh(r)\ 






ysh(r) 

















ch(r) / 



where ch(r) = cosh(2r) and sh(r) = sinh(2r), respec- 
tively. 

We observe that the nearest neighbor correlated GMPS 
requires only one squeezing parameter rg to generate the 
three-mode building block of Fig. [ljb). Furthermore, 
we can use finitely entangled TMS vacuum states with 
squeezing tt- For simplicity, we set rg = = r; n [23] 
and plot in Fig. [5] the squeezing strength needed to gen- 
erate the optimal input correlation <j)* ln for different noise 
correlations. For the Markov noise, in the plotted region 
the required correlation does not exceed 4>* n max ps 0.3, 
which can be realized by r injmax s» 1.08 (about 9.4 dB 
squeezing). For the non-Markovian noise, the required 
correlation does not exceed <j)* in max w 0.4, which cor- 
responds to T"in,max ~ 1-18 (about 10.2 dB squeezing). 
This shows that the required squeezing values for the 
presented setup could be realized with accessible non- 
linear media for a realistic assumption of noise correla- 
tions (these maximal squeezing values have recently been 
realized experimentally, see, e.g., Ref. |20|). 



IV. CONCLUSIONS 

We have demonstrated that a one-dimensional Gaus- 
sian matrix-product state, a multimode entangled state 
which can be prepared sequentially, can serve as a very 
good approximation to the optimal input state for encod- 
ing information into Gaussian bosonic memory channels. 
The fact that the GMPS can be prepared sequentially is 
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crucial because it makes the channel encoding feasible, 
progressively in time along with the subsequent uses of 
the channel. For the analyzed channels and noise mod- 
els, the GMPS achieves more than 99.9% of the Gaus- 
sian capacity and may be experimentally realizable as 
the required squeezing strengths are achievable within 
present technology. Furthermore, we have introduced 
a class of channel noises, originating from a chain of 
coupled harmonic oscillators at finite temperature, for 
which the GMPS is the exact optimal multimode input 
state. Given that GMPS are ground states of particu- 
lar quadratic Hamiltonians, our findings could serve as a 
starting point to find useful connections between quan- 



tum information theory and quantum statistical physics. 
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